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Finite-size scaling study of the d = 4 site-diluted Ising model* 
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We study the four dimensional site-diluted Ising model using finite-size scaling techniques. We explore the 
whole parameter space (density-coupling) in order to determine the Universality Class of the transition line. Our 
data are compatible with Mean Field behavior plus logarithmic corrections. 



1. Introduction 

If a pure system has a specific heat exponent 
a > the Universality Class changes when the 
dilution is introduced in the model (Harris crite- 
rion Jl|]), while it remains unchanged if a < (i.e. 
the Universality Class is that of the pure model) . 
In the Ising model in 4 (or 2) dimensions a = 
and the Harris criterion does not apply. 

Perturbative renormalization group (PRG) 
computations for d — 4 predict Mean Field with 
Logarithmic Corrections. On the other hand, pre- 
vious Monte Carlo (MC) results pointed to non 
Mean Field behavior In d = 2 there are also 
MC studies that conclude a change of the Univer- 
sality Class. 

We describe here the results of a higher statis- 
tics MC study @. A Finite-Size Scaling (FSS) 
approach has been used in order to study large 
lattices in the critical region. Results on the d = 2 
case are also briefly described j|. 

2. The model 

We work in a hypercubic four dimensional lat- 
tice. The action is: 



S = — ^2 ' '' i' T - a i- 



(1) 



<»j> 



where q are quenched uncorrelated random vari- 
ables whose value is 1 with probability p and 
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otherwise. 

For each {e^} configuration (sample) we per- 
form an Ising model simulation. 

There are two types of averaging. The first cor- 
responds to averaging in Ising configurations, and 
will be denoted with brackets, the second is as- 
sociated to the €i variables (sample average) and 
will be denoted by overlines. We first perform the 
Ising average, then the sample one. 

2.1. Observables 

For each spin configuration we measure the 
magnetization and first neighbor energy, defined 
respectively as 



M 



(2) 



We have focused our study in the following 
mean values (specific heat, susceptibility, Binder 
parameter, and correlation length, respectively) 
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where F is the Fourier transform of the magneti- 
zation at k = =£■. 
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3. Finite-size scaling techniques 

In a finite lattice at the the critical region, the 
FSS ansatz states that 



(0(L.J))=L^ 



o 



0(L~ 



(7) 



V L 

where lo is the corrections-to-scaling exponent, 
Fo is a (smooth) scaling function and xo is the 
critical exponent. For instance, x x = 7, x^ = u, 
and Xd p £ = v + 1. 

We study the quotient of 0{L\) and 0{L?) 



_ {Q(L 2 ,(3)) 

>0 = — — S 



Fo{^) 



-0(L-")(8) 



(0(Li,/3)) 

where s = j^. The unknown Fo can be elim- 
inated in this way: we look for j3{L\,L2) such 



that 



so 



+ 0(2T W ). 



(9) 



Now, this picture is slightly modified due to 
the presence of logarithmic corrections. Using the 
PRG analysis we obtain for the diluted model || 

£ oc L(logL)s 
dpi oc L 3 (logL)3 (i)V 2 V /I lF 



_ (10) 
X oc L 2 (log L) 4+tot 

C oc (logL^e-VSl 10 ^! 

The scaling behavior for the pure model is: 

£ oc L(logL)3 

X oc L 2 (logL)^ (11) 
C oc (logL)a 

It can be shown that £(£, (3,p)/ L remains as 
the scaling variable. 

The expected leading logarithmic correction 
for the v exponent goes as 1/ log L for the pure 
model, but changes to l/^/logX for the diluted 
case. For the r\ exponent, the correction is always 
order 1/ logL. 

4. Numerical methods 

We store individual measures of the energy and 
magnetization for extrapolating in a neighbor- 
hood of the simulation parameters. 
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Figure 1. Phase diagram of the d = 4 site-diluted 
Ising model. The tics are plotted at the sim- 
ulated values along the extrapolation direction. 
The black arrow indicates the percolation limit. 



The /3-derivatives are obtained from 

d W) = dp~W)=W)-WW), 



(12) 



which is a biased estimator. The statistical error 
behaves as 1/ y/Ns, N$ being the number of sam- 
ples and the bias as 1/Ni (Nj being the number 
of Ising independent measures in a sample). In 
our calculations \/N$ ~ Nj so a TVj — > 00 ex- 
trapolation is performed. 

As we gained statistics in a large number of 
samples with slightly different number of filled 
sites, in addition to a /3-extrapolation a p one is 
also possible. 

The probability of finding q site density for di- 
lution p is binomial. 

The observable value from a set of N$ samples 
at p' near p is 



<W,/3) = 
Ns fp' 
\P 



1-P' 
I-P 



It is also possible to compute p-derivatives but 
the statistical error is much larger than for the 
/3-derivatives (eight times typically). 

5. Results 

We use cluster methods to update the spin vari- 
ables. The diluted model is simulated in lattice 
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Table 1 

The v exponent for (L, 2L) pairs at the different 
dilutions. The last column corresponds to an in- 
verse linear extrapolation. 



V 


L = 8 


L = 12 


L = 16 


L = oo 


0.8 


.5175(11) 


.5154(11) 


.5142(13) 


.5110(25) 


0.65 


.5308(13) 


.5270(13) 


.5251(12) 


.5194(26) 


0.5 


.5482(16) 


.5428(19) 


.5412(19) 


.534(4) 


0.4 


.5604(15) 


.5532(19) 


.5478(18) 


.536(4) 


k, 0.3 


.5700(26) 


.5647(22) 


.5583(26) 


.549(5) 



the 



pure one 



sizes L < 2>2. We also studied 
(L < 64) as a check. 

We generate f 0,000 samples at p—0.8, 0.65, 0.5, 
0.4, 0.3. For each sample we take 100 nearly in- 
dependent measures after equilibration. 

Let first assume hyperscaling (there are not log- 
arithmic corrections). We use dpt; and \ t° obtain 
the critical exponents. In the v case we perform 
a L — > co extrapolation (see table @), using ui = 1 
(near the percolation value ||). For comparison, 
in the pure model we obtain v = 0.5019(14). 

For rj = 2 — j/i/ we find a weaker evolution with 
L and the dilution (see ref.||). 

We clearly discard the percolation values (y = 
0.686(2) and rj = -0.094(3) @) and a new fixed 
point for all the critical line is unlikely. 

Weak universality [y varying on the line) can- 
not be ruled out, but a more economic picture is 
a MF behavior plus logarithmic corrections (see 
% |). 

We can check the expected logarithmic correc- 
tions directly measuring £ and C at the critical 
point. To determine the critical coupling (or di- 
lution), we study the Binder 54 parameter. 

The hyperscaling violations for £ at the critical 
point are 



£(L,/3 c )cxL(logL)' 



(14) 



In table || we show the fitted 6^ values, which are 
reasonably close to the predicted values. 

6. The d = 2 model 

We perform 10 4 samples for L < 384 at dif- 
ferent dilutions. PRG predicts logarithmic cor- 
rections to the Ising behavior: v, equal to one, is 
corrected by 1/ logL and C oc log(log L). We find 
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Figure 2. Exponent v estimation for (L,2L) 
pairs. The solid lines are linear fits constrained 
to be v = 0.5 in the L — > 00 limit. 



Table 2 

5$ for the different dilutions. The second error 
bar is due to the uncertainty in the critical point 
and the first is the statistical one. 



V 


k 


P 


k 


1.0 


.198(3+5) 


0.8 


.203(5+10) 


0.65 


.181(6+14) 


0.5 


.189(7+12) 


0.4 


.241(8+12) 


0.3 


.281(10+23) 



also good agreement to the predicted behavior in 
both cases. 
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